Quantum Poincare Recurrences for Hydrogen Atom in a Microwave Field 
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We study the time dependence of the ionization probability of Rydberg atoms driven by a 
microwave field, both in classical and in quantum mechanics. The quantum survival probability 
follows the classical one up to the Heisenberg time and then decays algebraically as P(t) <x 1/t. 
This decay law derives from the exponentially long times required to escape from some region of 
the phase space, due to tunneling and localization effects. We also provide parameter values which 
should allow to observe such decay in laboratory experiments. 
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During the last two decades the manifestations of clas- 
sical chaos in microwave ionization of Rydberg atoms 
have been studied experimentally by different groups and 
many interesting results have been obtained In 
particular, laboratory experiments showed the quantum 
suppression of classically diffusive ionization process, in 
agreement with the predictions of dynamical localization 
theory [gj. The experimental technique based on accel- 
erated proton beams, which is used for production of hy- 
drogen atoms, allows to obtain interaction times with the 
microwave field of only few hundreds microwave periods 
On the contrary, the thermal beams used with 
alkali Rydberg atoms allow to vary the interaction time 
by orders of magnitude up to 10 5 microwave periods || . 
The first experiments of Walther's group || indicated an 
anomalously slow decay of the 10%-ionization threshold 
field as a function of the interaction time. This result 
cannot be explained within the picture of diffusive ion- 
ization in the domain of classical chaos. Some suggestions 
have been put forward to explain this slow decay which 
was attributed to some possible effects of noise for such 
long interaction times More recently, new experi- 

mental data for the behavior of the survival probability 
P(t) with time have been presented ||, showing an alge- 
braic law decay P(t) oc i~ Q , with a m 0.5. In the same 
paper, numerical simulations of quantum dynamics have 
been made, giving a value of a consistent with experi- 
mental data. The origin of the slow algebraic decay was 
attributed to the underlying structure of classical mixed 
phase space composed by integrable islands surrounded 
by chaotic components. However, the investigations of 
classical chaotic systems with mixed phase space showed 
that the probability of Poincare recurrences to the same 
region, or the survival probability up to time t, decays al- 
gebraically with power a w 1.5—3 [p|-pd|. Moreover, since 
the integral Jl P(r)dr is proportional to the measure of 



the finite chaotic region where the trajectory is trapped, 
the value of a should be greater than one. According 
to the correspondence principle, one expects that, in the 
semiclassical regime, classical and quantum systems ex- 
hibit the same decay law. Therefore the above exponent 
a w 0.5 found in the experiments requires, in this re- 
spect, an explanation. In particular the question arises 
whether this obtained value is generic or corresponds to 
some initial transient time behavior in a regime where 
quantum effects play an important role. 

Recent studies of quantum Poincare recurrences for the 
Chirikov standard map in the semiclassical regime with 
mixed phase space [|l2| showed that quantum P{t) fol- 
lows the classical decay during a relatively large time tu- 
The time tn gives the Heisenberg time scale, which is 
determined by inverse level spacings. For t > tji, the 
quantum survival probability starts to decay inversely 
proportional to time (a = 1) and becomes much larger 
than the classical one. The power a = 1 is due to ex- 
ponentially long times required to escape from some re- 
gion of phase space [Q. These exponentially long es- 
cape times are originated by tunneling from classically 
integrable region or by the exponential quantum local- 
ization. The above quantum behavior, with exponent 
a = 1, is different from the experimental data || and 
this constitutes an additional motivation for the present 
paper. Indeed the highly excited states of hydrogen atom 
in a microwave field can be described by the Kepler map 
which is very similar to the Chirikov standard map || 
and therefore one would expect the same behavior for 
the time dependence of survival probability. 

In order to investigate the probability decay for the hy- 
drogen atom in a microwave field, we choose the initial 
state with principal quantum number uq and numerically 
studied the survival probability P{t) in a linearly polar- 
ized monochromatic electric field e(t) — esin(wi). Here e 
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and u> are the strength and frequency of the microwave 
field, measured in atomic units. The quantum evolution 
is numerically simulated by the Kepler map ||, by the 
one-dimensional (lei) model of a hydrogen atom and by 
the 3d model for atoms initially prepared in states ex- 
tended along the field direction and with magnetic quan- 
tum number m = 0. The comparison of these three mod- 
els shows that the essential physics is captured by the Id 
model as already discussed in j^j. In addition we show 
that also the Kepler map gives an approximate correct 
description of the dynamics. 

The hydrogen atom in a linearly polarized monochro- 
matic electric field is described by the Hamiltonian 

p 2 1 

H=^- hezsin(wi), (1) 

2 r 

where, in the Id model the motion is assumed to take 
place along the field direction (z-axis, with z > 0). In 
order to compare classical and quantum dynamics it is 
convenient to use the scaled field strength eq = tn\ and 
frequency ujq = ujhq, which completely determine the 
classical dynamics. The classical limit corresponds to 
^off = ft/no —> 0, at constant cq, u>q. For u> > 1 the 
main change of the electron energy E occurs when the 
electron is close to the nucleus. As a consequence the 
dynamics is approximately given by the Kepler map |l. 

77 = N + k sin 0, 0= (j} + 2T:uj(-2ujN)- 3/2 , (2) 

where N = E/oj, k = 2.6ew -5 / 3 , <fi — tot is the phase of 
the microwave field when the electron passes through the 
perihelion and the bar marks the new values of variables. 
In the quantum case, the change of N gives the num- 
ber of absorbed photons while the number of photons 
required to ionize the atom is Ni — 1/(27JqOj). In clas- 
sical mechanics diffusive ionization takes place for fields 
above the chaos border: e > £c ~ l/(49wy 3 ) ||. The 
quantum dynamics of the model (|J) is described by the 
quantum Kepler map for the wave function ip(<f>): 

ip — exp(— iH (N))P exp(— ik cos fyip , (3) 

where H (N) = 2vr/V-2w7V, N = -id/ dtp, 4> = <t> 
(— oo < <fi < +oo), and the operator P projects prob- 
ability over the states with negative energy (N < 0) 
pL We introduce an absorption border for levels with 
n > n c , which for the Kepler map corresponds to N > 
N c « — l/(2nj?u>) O]. Such border occurs in real labo- 
ratory experiments, for example as a consequence of un- 
avoidable static electric field experienced by the Rydberg 
atoms during their interaction with the microwave field. 
The absorption border n c can be varied in a controlled 
way via a static electric field e s , the static field ionization 
border being e s n^, » 0.13. 

The results of quantum simulations for the situation 
similar to the experimental one (Fig. 2 (b) in Ref. ||) 



are shown in Fig. |lj. The Kepler map description allows 
us to study the quantum dynamics up to very long times 
(t = 10 8 , here and below time is given in microwave 
periods). In the case no = 23 the quantum data for 
the survival probability P(t) obtained from the quan- 
tum Kepler map and the l<i hydrogen atom model agree 
with each other (see the inset of Fig. [l]) and with the 
numerical computations of ||. However all these data 
are strongly different from the classical probability decay 
shown in Fig.|l], which displays a slope a w 2. The rea- 
son of this disagreement should be attributed to the fact 
that no — 23 is not in the semiclassical regime. Our data 
for the Husimi distribution, obtained from the Wigner 
function by smoothing over the size Ti | fl4| , show that a 
significant part of the probability is trapped inside the 
stable island at n « 20 (un 3 w 1). For this reason the 
probability decays slowly during a long time t w 10 5 after 
which it drops faster. If no is increased significantly, the 
semiclassical regime is reached and the quantum proba- 
bility decay becomes close to the classical one up to the 
time scale tn ~ 10 . Our data show that tn is propor- 
tional to no (at fixed eo> w o)j m agreement with previous 
estimates of Ref. p2[ , according to which tn oc l/h e s. 
After this time the quantum 1/t decay is clearly observed 
in agreement with the results of ]l2|. 
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FIG. 1. Survival probability P(t) as a function of the inter- 
action time t (in units of microwave periods) for eo = 0.065, 
wo = 1.6, n c — 2no: quantum Kepler map for no = 23 (thin 
solid line), no = 230 (thin dotted line), no — 2300 (thick 
solid line) and classical Kepler map (thick dashed line, en- 
semble of 10 9 trajectories). The straight lines have slopes 1 
and 1.94, the latter coming from a fit of the classical decay 
for 2 x 10 2 < t < 3 x 10 4 . Inset: quantum Kepler map (solid 
line) versus Id hydrogen model (dot-dashed line) for no = 23. 

In Fig. we show a more realistic case in which, ini- 
tially, classical and quantum probabilities decay in a very 
similar way and where only after a time tn ~ 5 x 10 2 , 
the quantum survival probability starts to decay more 
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slowly (P(t) oc 1/t) than the classical one which decays 
as l/t a , with a « 2.15. This case corresponds to no = 60 
and can be observed in experiments similar to those per- 
formed in S. Again the quantum Kepler map gives a 
qualitatively correct description of the ionization process 
up to very long interaction times. The comparison of 
quantum simulations for the Id hydrogen atom model 
and the 3d dynamics is shown in the inset of Fig. |^. It 
demonstrates that both dynamics give very close results, 
confirming that the essential physics is captured by the 
Id model. We put the absorption border near the initial 
state (n c = 64) in order to have p c = £<t,/AN c w 3.5 > 1, 

where 1$ = 3.3e§w^ 10 ^ 3 no is the localization length in 
number of photons || and AN C = {no/2u)o)(n^/n 2 c — 1) 
is the number of photons required to reach the absorp- 
tion border. In this way the probability can go out very 
easily and the 1/t probability decay is observed after a 
short transient time of the order of 20 microwave peri- 
ods. On the contrary, when p c < 1, as in the case of Fig. 
|l| for no = 23 (p c « 0.3), strong fluctuactions around 
the 1/t decay take place. This is analogous to the huge 
(log-normally distributed) conductance fluctuactions in 
a disorder solid with localization length smaller than the 
system size [|l5| . 
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FIG. 2. Survival probability for eo = 0.1, uo = 2.6, 
no = 60, n c — 64: quantum Kepler map (dotted line), 
quantum (solid line) and classical (dashed line, ensemble of 
3 x 10 6 trajectories) lei hydrogen model. The straight lines 
have slopes 1 and 2.15, the latter coming from a fit of the 
classical decay for 5 x 10 2 < t < 2 x 10 4 . Inset: quantum 
survival probability for the Id model (solid line) and the 3d 
model (dot-dashed line). 

In order to confirm that the algebraic probability de- 
cay is related to the sticking of classical trajectories and 
of quantum probability near the integrable islands in the 
phase space, we show in Fig. |^ the time evolution of 
the survival probability distribution in the phase space 
of action-angle variables (n, 8) for the Id model. In the 



classical case 3 x 10 6 orbits were initially homogenously 
distributed in the angle 6 on the line n = hq = 60, cor- 
responding to the initial quantum state with principal 
quantum number n = 60. After 50 microwave peri- 
ods, the classical distribution of non ionized orbits shows 
a fractal structure which surrounds the stability islands 
(Fig. U top left). At larger times this distribution ap- 
proaches more and more closely the boundary critical in- 
variant curves (Fig. || bottom left). One of them confines 
the motion in the region with n > rib ~ "o(ec/eo) 1 ^ 5 ~ 41 
where rib determines the classical chaos border for given 
eo- Other invariant curves mark the critical boundaries 
around internal stability islands (for example at n ~ 55, 
corresponding to tun 3 rj 2). In the quantum case the 
value of h e g is not sufficiently small to resolve the fractal 
structure at small scales. However, the Husimi function 
shows similarities with classical probability distribution 
at t = 50 (Fig. H top right). At longer times, the diffu- 
sion towards the boundary at rib is slowed down due to 
localization effects and penetration of the quantum prob- 
ability inside the classical integrable islands. At t = 10 4 
(Fig. U bottom right) the quantum probability is concen- 
trated in a layer near rib- Due to localization effects, the 
Husimi function does not change significantly for a very 
long interaction time (10 3 < t < 3 x 10 4 ). Eventually 
the probability starts to penetrate very slowly inside the 
main island at n ~ nj. Therefore tunneling and localiza- 
tion effects are responsible for the slow 1/t decay of the 
quantum survival probability seen in Fig. ^. 

The fractal structure of the classical distribution is 
washed out at scales smaller than the minimal quantum 
cell h Q g. Therefore a better resolution can be obtained 
increasing the principal quantum number no, at fixed eo, 
ujq, and n c /riQ. The Husimi function clearly reflects the 
underline fractal structure at very high principal quan- 
tum numbers n = 150 (Fig. |J left) and n a = 1200 (Fig. 
|] right). Similar quantum fractals have been found in the 
kicked rotator model with absorbing boundary conditions 

s 

Notice that the probability decay P(t) is related to 
correlations decay via C(t) oc tP(t) ||. In the case 
of a = 1 this implies that correlations do not decay. 
The Fourier transform of C(t) gives the spectral density 
S(lj) of the effective noise produced by the dynamics: 
S(u>) — J C(t) exp(iu>t)dt « 1/ui. This shows that the 
spectral noise associated with the quantum Poincare re- 
currences with a = 1 scales like S(uS) oc 1/uo. A similar 
behavior of noise has been observed in many scientific dis- 
ciplines Jit] ], for example in the resistence fluctuactions 
of different solid state devices Es]. This phenomenon 
is known as 1/f noise and usually extends over several 
orders of magnitude in frequency, indicating a broad dis- 
tribution of time scales in the system. In the case of 
quantum Poincare recurrences this property stems from 
the exponentially low escape rate from some regions of 
the phase space. 
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In summary, on the basis of our previous investigations 
and of the numerical studies presented in this paper we 
conclude that the survival probability for Rydberg atoms 
in a microwave field decays, up to the time scale tn oc no, 
in a way similar to the classical probability. For t > tji 
the quantum probability starts to decay slower than the 
classical one, with the exponent of the algebraic decay 
a = 1. We have given parameter values which should 
allow to observe quantum Poincare recurrences in mi- 
crowave experiments with Rydberg atoms. 




FIG. 3. Classical density plot (left) and Husimi function 
(right) in action-angle variables (n, 0), with 30 < n < 63.5 
(vertical axis) and < 9 < 2tt (horizontal axis), for the Id 
model in the case of Fig. 2. Husimi function is averaged in 
a finite time interval to decrease fluctuations: 50 < t < 60 
(top); 2 x I0 3 < t < I0 4 (bottom left); 9.9 x I0 3 < t < I0 4 
(bottom right). The color is proportional to the density: blue 
for zero and red for maximal density. 




FIG. 4. Quantum fractal Husimi function for parameter 
values and interaction times as in Fig. ^ top left, with 
no = 150 (left) and n = 1200 (right). 

This research is done in the frame of EC program 
RTN1-1999-00400. 
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